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Abstract 

We consider a dynamical system subjected to weak but adiabatically slow 
fluctuations of external origin. Based on the "adiabatic following" approximation 
we carry out an expansion in a|/x| _1 , where a is the strength of fluctuations and 
l/il -1 refers to the time scale of evolution of the unperturbed system to obtain a 
linear differential equation for the average solution. The theory is applied to the 
problems of a damped harmonic oscillator and diffusion in a turbulent fluid. The 
result is the realization of 'renormalised' diffusion constant or damping constant 
for the respective problems. The applicability of the method has been critically 
analyzed. 
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I. Introduction 



The standard paradigm of the temporal evolution of nonequilibrium processes re- 
garded, in general, as stochastic processes is the century-old problem of Brownian 
motion [1,2]. This involves the random motion of microscopic particles effectively in- 
troducing the motion of a physical system, the Brownian particle to be observed on a 
macroscopic level. To generate the successive levels of description from the microscopic 
to the macroscopic realm one essentially introduces coarse-graining of space and time 
in the dynamics. Although there exists no general program of coarse-graining, it is 
nevertheless possible to realize the dynamics of stochastic processes in terms of some 
systematic separation of time scales consistent with the experiments at the macroscopic 
level of description. 

The standard separation of time scales in the description of Brownian motion in- 
volves correlation functions which are nonzero over some interval r c which is the cor- 
relation time of fluctuations and we require that At, the coarse-grained timescale over 
which one observes the average motion is much greater than r c , such that 7 _1 ^> At 3> 
r c , where 7 _1 is the system's damping time. Physically this implies that one smoothes 
out the fluctuations of the system on a time scale during which microscopic particles 
are correlated but not on a scale during which the system is damped. Thus the fluc- 
tuations considered in the stochastic process of the Brownian motion are weak and 
rapid. 

In the present problem we consider a multivariate dynamical system driven by weak 
but adiabatically slow fluctuations. The slow fluctuations characterized by very long 
correlation time have also attracted a lot of attention of various workers over the years 
[2-4,7,8,11,14,15]. While the overwhelming majority of the treatment of stochastic 
differential equation with fast fluctuations is based on the assumption that there is a 
very short auto-correlation time r c , such that one can adopt the scheme of expansion 
in ar c , suitable simplifying approximation for dealing with very long correlation-time 
is relatively scarce. In general, the problem of long correlation time is theoretically 
handled at the expense of severe restriction on the type of stochastic behavior. For 
instance, several authors [2-4,7,8,11,14,15] have tried the linear and nonlinear models 
within the framework of Markov processes of the type dichotomic processes, two-state 
Markov processes, random telegraphic processes etc. Our aim here is to explore a 
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perturvative method for finding an equation for the average solution pertaining to the 
separation of timescale implied in the inequality 

1 

T-r < At < T c , 
W 

|// 1 being the largest eigenvalue of the unperturbed system, where we do not keep any 
restriction on the type of stochastic behavior. The strategy being perturbative is based 
on an expansion in al/il" 1 rather than in ar c as is done in the case of fast fluctuations. 
The method dealt in the present treatment is thus somewhat complementary to the 
scheme of expansion of the latter kind. 

To put the issue in a proper perspective we first borrow a simple example of adiabatic 
dynamics in terms of Bloch equations [5,6], wellknown in magnetic resonance and 
quantum optical experiments. The problem concerns a two-level system interacting 
with a single mode electromagnetic field, where the field £{t) varies slowly enough 
"adiabatically" in the time scale of the inverse of the damping constant or frequency 
detuning between the atom and the field. The term "adiabatic following" is thus used 
to describe collectively the associated experimental phenomena [5,6]. The model is 
described by the following equation, 



/ 


u 


\ 


( 


1 

T 2 


-A 





\ 


/ 


U 


\ 




( 





\ 




V 






A 


_ j_ 


gS(t) 






V 




+ 









V 


w 


J 


V 







i 




V 


w 


J 




V 


W c<1 

Ti 


/ 



Here u, v, w are the Bloch vector components, T\ and T 2 are the energy and dephas- 
ing relaxation terms, A is the detuning of the frequency of the field £{t) from that of 
the two-level system, g includes the effect of coupling of the atom to the field. The 
equilibrium value, towards which the population inversion w relaxes when S = is 
denoted by w eq . Adiabatic following approximation asserts that if the field S(t) is var- 
ied slowly enough then w, the population inversion variable would follow adiabatically 
from —1 to ~ +1 in the process, i.e., a ground state population could be adiabatically 
inverted. 

Our problem in the present investigation concerns such processes where the adiabatic 
variation of £(t), in addition, is stochastic. Thus the usual limit in the "adiabatic 
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following" applies, i.e., the rate of variation of the pulse or fluctuations is much small 
compared to the relaxation rate of the system. With these in mind we may treat 
Eq.(l) as a stochastic differential equation provided the stochastic properties of £(t) 
are a priori known. 

To formulate the problem we thus consider a system subject to fluctuating external 
forces where the fluctuations are weak and adiabatically slow. The equation of motion 
then become a stochastic differential equation, a particular category of which is a 
general form of Eq.(l) [For simplicity we disregard the constant part on the right hand 
side] , 

ii = A(t)u. (2) 

Here A(£) is a random function of time, stochastic properties of which are given. 
Linear multiplicative noise ( Eq.(2) ) has got wide application in studying the random 
Markov process [7], fluctuating barrier crossing [8], enzymatic kinetics in biology [9], 
nuclear magnetic resonance in physics [10] and stochastic resonance in linear system 
[11] and in many other context [12]. 

Based on the systematic separation of time scales using adiabatic following approx- 
imation a differential equation for the average solution (u) is obtained. This approx- 
imation allows us an expansion in where a is a measure of the strength of 
fluctuations and /z -1 refers to the internal time scale of the unperturbed system. 

As an immediate application of the method we treat the problem of a damped 
harmonic oscillator with adiabatically fluctuating frequency. The method is extended 
to the problem of diffusion in a turbulent fluid as another illustration. The central 
result is that one realizes a "renormalised" transport coefficient or a damping constant 
so that the diffusion or the damping process gets significantly modified by adiabatic 
stochasticity. We show that the method is equipped to deal with similar kinds of 
stochastic processes. 

The outlay of the paper is as follows; In Sec. II we discuss the method of adiabatic 
following approximation on the stochastic differential equation of the form (2). The 
essential idea is to extract the average dynamics of the relevant physical quantity. The 
method has been critically analyzed in Sec. III. The method is applied to two specific 
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cases in Sec. IV. We point out that a wide class of problems can be treated in a similar 
way. The paper is concluded in Sec.V. 

II. A method for weak and adiabatic fluctuations 

To start with we consider a linear equation of the type (2) and rewrite it as 

ii = { A + aAi(t)}u , (3) 

where u is a vector with n components, A is a constant matrix of dimension n x n 
and Ai(£) is a stochastic matrix, a is a parameter ( of dimension 1/t) which measures 
the strength of fluctuation. 

It is convenient to assume that Ai(i) is a stationary process with (Ai(£)) = 0. 
Eq.(3) admits of two important time scales of the system measured by the inverse of 
the largest eigenvalue of the matrix Ao and the time scale of fluctuations of Ai(t) 
(the correlation time of fluctuation). As already been mentioned that in the treatment 
of overwhelming majority of the stochastic processes, such as, motion of a Brownian 
particle in a fluid or electromagnetic waves in a turbulent atmosphere, one essentially 
considers a situation where the fluctuations are weak and rapid. The correlation time 
of fluctuations is much short compared to the time scale set by the inverse of the 
eigenvalues of A . In the appropriate limit we encounter the delta-correlated events 
and solve approximately or exactly the relevant stochastic differential equations [2]. 
The familiar examples of paramagnetic resonance and line broadening are well known 
in this context. 

Since in the present problem we consider a stochastic process in which the fluctua- 
tions are weak but adiabatically slow, Ai(£) is an adiabatic stochastic process. There- 
fore the usual procedure of systematic cumulant expansion which inherently takes into 
account of short correlation time of fluctuations is not valid. An alternative treatment 
is thus sought for. 

To this end we first introduce an interaction representation as given by, 
u(t) = exp(A t)v(t) 
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and applying it to Eq.(3) we obtain, 

v = aV(t)v , 
where, V(t) — exp(— A t) Ai(t) exp(A £) . 

On integration the last equation yields, 

v (t) = v (0) +a f V(t')v(t')dt' . (4) 
Jo 

On iterating the Eq.(4) once, we are now led to an ensemble average equation of the 
following form, 

( v (t)) = v(0) + a 2 f dt' C dt"(V(t')V(t")v(t")) . (5) 
JO Jo 

The equation is still exact since no second order approximation (as usually done) 
has been used. 

Now taking the time derivative of v(t) we arrive at the following integrodifferential 
equation in which the initial value v(0) no longer appears, 

j t (v(t)) = a 2 j\v(t)V(t')v(t'))dt'. (6) 
Making use of a change of integration variable i! = t — r we obtain, 

d ( v (t)) = a 2 [\v(t)V(t - r)v(t - r))dr . (7) 



dt Jo 
Reverting back to the original representation Eq.(7) yields 

d 



— {u{t))=A (u) + a 2 f (Ai(i)exp(Aor)Ai(t-r)u(t-r))dr. (8) 
dt Jo 

The adiabatic following assumption, that Ai(t) and the components of u(t) vary 
slowly on the scale of inverse of A , can now be utilized. Following Crisp [6] we note 
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that a Taylor series expansion of Ai(t — r)u(t — r) in the average (. . .) of the a 2 -term in 
Eq. (8) allows the integral to be evaluated and the last equation reduces to the following 
form, 



Jt Ht)) = A (u) + a 2 Y: o [ -^r(Mt){l rfrexp(A r)r"}— [A^u®]) .(9) 

The integral in Eq.(9) can be evaluated by rewriting it in terms of the following 
matrix elements 



ft 

EA 3 ^p^\ Re //,, < 



Therefore, 



T = n!DE n+1 D" 1 , 



(10) 



where D is a matrix which diagonalises A and 
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and jJL jj £1X6 the eigenvalues of A . 

The Eq.(9) then assumes the following form, 



j t (u(t)) = A Q (u) + a^i-lTiA^BE^B-^iMitHt)]) . (11) 
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Although the equation involves an infinite series it is expected to yield useful result 
in the adiabatic following approximation. If this approximation is valid, the quantity 
[Ai(i)«(i)] varies little in the time l/^™/ 1 ! -1 of E n+i and the series converges rapidly. 
Keeping only the two lowest order terms we arrive at, 

^-(u(t)) = Ao(«) + o^A^XxA^Mt)) - a 2 (A 1 (t)X 2 A 1 (t) W (i)) 
at 

-a 2 (A 1 (t)X 2 A 1 (t)u(t)) (12) 

where, 



X^^DE^D- 1 . (13) 

It is evident that the average (u) is related to a more complicated average. As a 
next approximation [13-15] we now suppose that the latter averages may be broken up 

as, 



(A 1 (t)X 1 A 1 (t)u(t)) » (A 1 (t)X 1 A 1 (t))(u(t)) (14) 
and so on. Keeping only the terms upto the order of a 2 we obtain, 



-(u(t)) = {a + ^[(Ax^XxAx^)) - (A 1 (t)X 2 A 1 (t)) 

-{A 1 (t)X 2 A 1 (t))A ]}{u(t)). (15) 

Thus the average of u(t) obeys a nonstochastic differential equation in which the 
effect of weak adiabatic fluctuations is accounted for by "renormalizing" A through 
the addition of constant terms of the order of a 2 . The net effect is that depending on 
the specificity of the situations one realizes a dissipative or a gain term in the average 
dynamics. We note in passing, that the average dynamics of u is independent of any 
explicit correlation function. 



III. Discussions on the method 
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Theory of stochastic differential equations with multiplicative noise has a long his- 
tory. The stochastic processes dealt with the overwhelming majority of the cases con- 
cern the fast processes (more precisely, the correlation time between the noise events 
has been considered to be the shortest timescale of the dynamics). In the previous sec- 
tion we have considered a stochastic process which adiabatically slow. The traditional 
scheme of solving stochastic differential equations with fast noise processes is that one 
reduces them to Bourret's integral equations [13] and then performs the decoupling of 
the product of operators. In the present paper we have followed equation scheme upto 
Eq.(9) and then make use of "adiabatic following" approximation. It is necessary to 
make the following distinctions; 

First, note that in going from Eq.(8) to (11) we have made no approximation so far 
as the full infinite series is concerned. Also each term is not of order ar c as in the case 
of fast processes (as emphasized by van Kampen [14]) but of order a J^-[A(t)-u(t)]//i^ +1 . 
Just as the theory of fast processes is valid for ar c very small which implies that the 
successive cumulants in the expansion are small, validity of the description of adiabatic 
processes rests on the smallness of successive ajj^[A(t)u(t)]/ fJ>]j~ l terms. Thus the two 
expansions are essentially different. 

Second, the decoupling approximation has been carried out both in the fast as well as 
in the slow processes. Its justifications in the former case has been established early by 
Brissaud and Frisch [15]. It had been strongly advocated by van Kampen [14] who has 
asserted that although it seems to neglect certain correlations, the "statistical mechan- 
ics of transport processes would be in a very sorry state without such hypothesis" . It is 
not difficult to comprehend that its spiritual root lies in "stosszahlansatz" , "molecular 
chaos assumption" or "random phase approximations". The essential point, however, 
in the decoupling scheme is the realization of a separation of time scale of average of 
the product of fluctuating quantities A(t) and the average of u itself consistent with 
the expansions pertaining to slow or fast processes. 

That the two expansion schemes in the fast and slow stochastic processes are different 
can be confirmed if one compares the lowest order terms of the corresponding evolution. 
According to the present scheme Eq.(15) itself asserts that (free motion neglected) 
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where 1 refers to the timescale set by the A matrix which is short in the adiabatic 
following limit. For a fast stochastic process on the other hand the counterpart of 
Eq.(16) is 

j t (u)~a\(u) , (17) 

where r c defines the very short correlation time of the noise [14]. 

The difference in the expansion schemes also makes the relative errors made in the 
decoupling approximation in the two cases, different. To this end we first note that 
Eq.(15) is obtained from Eq.(8). Upto second order it means omitting terms of the 
order (a At) 3 and higher (where At defines the coarse-grained time scale of evolution 
of the average). As the lower bound of At is determined by it implies that we 

neglect terms of the order (a|yu| -1 ) 3 in the evolution equation. Thus the relative error 
made in the decoupling approximation in the case of adiabatic expansion is (a|/i| -1 ) 3 . 

As demonstrated by van Kampen [14] the corresponding error made in the decou- 
pling approximation in Bourret's scheme is of the order (ar c ) 3 . The workability of the 
decoupling approximation in the fast and slow stochastic processes is thus demonstrated 
in the two different expansion procedures ensuring their respective fast convergence in 
the limit ar c (fast processes) or a^l -1 (slow processes) small but finite. 

So, to summarize, we point out that the implementation of Bourret's decoupling 
approximation is a major step for almost any treatment of multiplicative noise upto 
date [2,7,11-15]. This is because of the fact that the average of a product of stochastic 
quantities does not factorize into the product of averages, although it has been argued 
that [2,7,11-21] good approximations can be derived by assuming such factorization. In 
the case of fast fluctuations it has been justified if the driving stochastic noise has a fast 
correlation time such that Kubo number a 2 r c is very small in the cumulant expansion 
scheme ( an expansion in ar c ). The factorization has been shown to be exact in the 
limit of zero correlation time and in some specific noise processes [7,14] and the solution 
for the average can be found exactly. The present scheme of adiabatic expansion on 
the other hand is an expansion in a^l -1 and it may be argued in the same way that 
factorization in the slow fluctuation is valid where o; 2 |/x | 1 is very small. Essentially 
it implies that u(t) in the average ( in the right hand side of Eq.(12) ) is realized as 
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an average (u(t)) ( which varies in the coarse-grained timescale At ) pertaining to the 
separation of the timescales in the inequality \fi\~ 1 C At < r c adopted in the present 
case instead of the inequality r c C At C employed in the case of fast fluctuation 
and cumulant expansion. 

IV. Applications 

A. Damped harmonic oscillator with adiabatically fluctuating frequency 

To illustrate the above-mentioned method we consider a model of damped harmonic 
oscillator with random frequency where the fluctuation is weak and adiabatically slow 
in the time scale of dissipation. The opposite limit of weak and rapid fluctuations in 
frequency has been studied by numerous authors in connection with turbulence, wave 
propagation, line-broadening [10], lasers, chaotic dynamics [20,21]. A comprehensive 
treatment has been given in van Kampen [14]. 

We are now in a position to apply the result(15) to the following equation, 

x + u 2 (t)x + -fx = , (18) 
with an adiabatically stochastic frequency, 

cu 2 (t)=u 2 [l+am, (19) 

where is an adiabatic stochastic process with zero mean (£(£)) = 0; c^o is the 
frequency of the unperturbed oscillator and 7 is the damping constant, a, the smallness 
parameter is dimensionless in Eq.(19). 

Rewriting Eq.(18) in the form, 
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one identifies, 



-u)q -7 



1 \ 9>/ . / \ I x 



(21) 



of Eq.(15). Xx and X 2 are related to E x and E 2 through (13) and are given by, 
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where A and B are related to the eigenvalues (ei, e 2 ) of A matrix : 



A = ~l + y^-Au;l ( ei ), 
5 = | + ^\/7 2 -^o 2 , (-e 2 ). 



D matrix is given by , 
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A -B 
\ (A+l)2 (B+l)i / 



So for the present problem, the Eq.(15) takes the form, 
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d ( (x)\ ( 1 \ / (x) \ 

dt [ (y) ) = [ -^-a^ (e)-a^O - 7 + « 2 ^ 2 )7 ) { (v) ) ' 

or 



(x) + 7 [1 - a 2 (e)](x) + K 2 + cA; 2 <£ 2 > + a^mix) =0. (22) 

It is thus evident that the adiabatic fluctuations in frequency cause a suppression of 
damping of the average amplitude of the oscillator. Or, in other words, the dissipative 
oscillator experiences a partial gain in average amplitude by an amount, 



7gai„ = a 2 7 (£ 2 W>. (23) 

As expected, the frequency of the unperturbed oscillator has also undergone a shift 
in addition to this gain in amplitude. 

The above result can be compared to the case of fast fluctuations in frequency as 
treated by van Kampen and others. It is interesting to note that while the adiabatic 
fluctuations result in gain in amplitude the fast fluctuations cause a damping of the 
average amplitude, in general. This damping may even be negative when the fluctua- 
tions are particularly strong at twice the unperturbed frequency. The latter results had 
been found to be useful in the context of a fluctuation-dissipation relation in chaotic 
dynamics [21]. 

The theory developed in the preceding section also permits us to calculate the dy- 
namics of the higher moments. For example, the equations of the three moments can 
be found from Eq.(18), 
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or rewriting it in the form, 
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we identify, 
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The eigenvalues of A are, 



and Ai(t) = aw 2 £(t) 
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(26) 



e 3 = -7 - V7 2 - 4 ^o 
Eq.(15) therefore takes the form of an evolution equation of higher moments, 
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where c\ = (£ 2 ) and c 2 = . 

What follows as a consequence of Eq.(27) is the shifting of eigenvalues of the un- 
perturbed oscillator. The eigenvalue which corresponds to e\ of the unperturbed case 
now becomes, 

2 J 7 1 \ 

-7 - a < -ci + — -2 ? ci - -c 2 > 

to second order in a. Thus the damping of energy of the unperturbed oscillator gets 
enhanced beyond a critical value depending on the positivity of the term included in 
the parenthesis of the last expression. In the negative region the term acts as a gain 
term leading to a suppression of dissipation of energy of the oscillator. 

In contrast to this case of adiabatic fluctuations, fast fluctuations make the oscillator 
unstable energy-wise due to the fluctuations in the force that have twice the frequency 
of the oscillator. Addition of damping term, however, may result in stability below a 
certain critical value and unstability above it. This result has been particularly relevant 
in establishing the Kubo relation in chaotic dynamics [20]. 

B. Diffusion in a turbulent fluid 

As a next application we consider a diffusive process in a fluid in motion described 

by, 

^ = -V-H + W 2 n. (29) 
at 

Here n(r,t) is the number of 'probe particles' per unit volume and v(r, t) is the 
velocity of the moving fluid. If turbulence sets in then v(r, t) becomes a stochastic 
function of r and t. The problem has been addressed by numerous workers over several 
decades. To quote a representative few of them we refer to Ref [14]. 

For the present problem we consider v(r,t) as an adiabatic stochastic process. The 
problem then is to find an average n(r, t) for the given initial condition 
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n(r,0) = S(r). 

We consider the turbulence to be weak and slow and without any loss of generality 
assume 



(v(r,f)> = 

Eq.(29) is of the form (3), provided the matrix A and Ai correspond to 

A = DV 2 and aAi = V • v 

The symbol V, as usual, acts on every functions of r that appears to the right of it. 
Eq.(8) then takes the following form, 

^(n(r, t)) = DV 2 (n(r, t)) + jT dr(V • v(r, t)e^V ■ v(r, * - r)n(r, t - r)> . (30) 
We take the Fourier transform in space of the last Eq.(30) to obtain, 



-(n(k,t)) = -Dk 2 (n(k,t))-(2n)~^ dr / ( W< (q, t) Vj -(q', t - r) 



n(k,t - r)> x (ft + ^ + fcije-^^') 2 ^. + fc^qdq'. (31) 

Expanding i>j(q', i — 7")n(k, i — r) as a Taylor series and integrating over r as before 
we arrive at, 



£<n(k, t)> = -Dk 2 (n(k, t)) - (2vr)- 3 £ £(-!)" / / <T*(q, 

OT ■ • „ J -co J -co 



i,j n 

1 ^>,(q',t)n(k,t)]> X (ft + ft' + kM + %)rfqrfq'. (32) 



[D(k + q') 2 ] n+1 
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Taking into account of the following properties of the stochastic function v(q, t), 



(^(q,0^(q , ,i)) = 5(q+q')(27r)iC ij (q), 

(^(q,%(q',t)> = 5(q + q')(27r)f^.(q), (33) 

and imposing the adiabatic following approximation that Vj(q, t — r)n(k, t — r) varies 
much slowly in the time scale of D~ l we obtain, 

f) { , roo i 

-(n(k,t)) = -DA; 2 - (2tt)-> £ j^Q^Hkj - g,-) - (k )2 <fr 

o ^ /'OC I 

+(2t)" 5 E /_ ^(q)^(% - ^) D2(k _ q)4 ^q 

o roo i-2 ! 

-(27T)-' E ^(q)^(% - D{k _ q)4 ^q| • (34) 

The effect of adiabatic stochasticity in the motion of the fluid thus essentially is 
to recover a renormalised diffusion coefficient D(k) of 'test' particles in the following 
form, 

3 f°° 1 

D(k) = d + E ^(q)fci(fc 3 - - gj) Dfc2(k _ q)2 ^q 

o /-oo I 

-(*)-' £ /_„ c'MMk, - 

3 /'OO ] 

+(2ir)~J2 C^qMkj - qj )— -dq. (35) 

J-oo Lf{k - qj 

Hence Eq.(34) reduces to, 

^(n(k,t)) = - J D(A;)A; 2 (n(k,t)) (36) 

a normalized diffusion equation for isotropic turbulence. In the absence of any detailed 
knowledge about the stochastic properties embedded in CV,-(q) and C^-(q) it is diffi- 
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cult to proceed further. Nevertheless in the limit of small q one might expect some 
interesting behavior as has been observed in the case of rapid fluctuations. 

It is well known that the long wavelength fast fluctuations are insufficiently damped 
by the viscosity (which appears as a parameter in the correlation function of the 
incompressible fluids) which ensures the existence of a finite r c . This causes long 
time tails in the correlation functions. As van Kampen [14] emphasized the stochastic 
description in terms of an average (n) ceases to become meaningful in these cases. Since 
the present treatment of slow fluctuations is free from explicit correlation functions and 
we deal only with averages, such pathological problems of long time tails or memory 
need not trouble us to that extent. This leads us to believe that the average description 
remains more meaningful in such cases. 



V. Conclusion 



In this paper we develop a method for treatment of weak but adiabatically slow 
stochastic process. Based on the "adiabatic following" approximation we recast a class 
of linear stochastic differential equations with multiplicative noise into a differential 
equation for the average solution. This has been carried on the basis of an expansion 
in where a is the size of the fluctuation and refers to the time scale of 

evolution of the unperturbed system. The result differs significantly from the corre- 
sponding treatment of weak and rapid fluctuations which relies on the expansion in 
ar c , where r c is the auto-correlation time of fluctuations [14]. It is also necessary to 
emphasize that in the present work no a priori assumption on the nature of noise in 
Ax(t) (like Ax(t) is a Gaussian random process etc., which has received so much atten- 
tion in the literature) has been made. Although in our applications we have dealt with 
classical and linear problems, the theory can be extended to quantum mechanical and 
nonlinear problems as well. We hope to address such issues in future communications. 
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